Estimates of Effective Hubbard Model Parameters for C20 isomers 
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We report on an effective Hubbard Hamiltonian approach for the study of electronic correlations 
in C20 isomers, cage, bowl and ring, with quantum Monte Carlo and exact diagonalization methods. 
The tight-binding hopping parameter, t, in the effective Hamiltonian is determined by a fit to 
density functional theory calculations, and the on-site Coulomb interaction, U /t, is determined by 
calculating the isomers' affinity energies, which are compared to experimental values. For the C20 
fullerene cage we estimate fcagc ^ 0.68 — 1.36 eV and (f7/t) ca ge — 7.1 — 12.2. The resulting effective 
Hamiltonian is then used to study the shift of spectral peaks in the density of states of neutral and 
one-electron-doped C20 isomers. Energy gaps are also extracted for possible future comparison with 
experiments. 

PACS numbers: 71.20.Tx, 71.10.Fd, 02.70.Uu 



I. INTRODUCTION 

The successful synthesis of gas phase C20 molecules 
displaying the dodecahedral fullerene cage structure^ has 
induced considerable research interest in this smallest 
fullerene molecule, partly because of the superconducting 
property of electron-doped Cgo materials^ in the same 
fullerene family, and partly because previous theoreti- 
cal speculations^^ about the existence of such a cage 
molecule were confirmed. In the experiment three ma- 
jor C20 isomers, of cage, bowl and ring structures, were 
produced. See Fig. Q]for the molecular structures. Pho- 
toelectron spectra (PES) were also measured for all three 
brominated and one-electron doped isomers (C^). Affin- 
ity energies (AE, see definition below) were then ex- 
tracted from the PES figure, giving ^4_E cagc = 2.25 eV, 
AE how \ = 2.17 eV, and AE ring = 2.44 eV, respectively. 
Since PES spectra reflect both the isomers geometric 
character and the strength of electronic correlation inside 
the molecules, it is a unique opportunity to investigate 
the interplay between geometry and electronic correla- 
tions in the three isomers. Here, we shall do this by es- 
timating the parameters in an effective Hubbard model 
description of these isomers, parameterized in terms of 
the on-site repulsion, U, and hopping integral, t. We 
find that for the geometry with the highest curvature, 
the fullerene cage, correlations effects as measured by U/t 
are relatively important. Experiments have also shown 
evidence for solid phases of C20 fullerene cageS^ further 
emphasizing the need for understanding strong correla- 
tion effects in this isomer. 

The possibility of superconductivity arising from a 
purely electronic mechanism in the Cgo fullerenes was 
suggested in the early 90's£ and supported by pertur- 
bative calculations^ starting from a one-band Hubbard 
model. In this picture superconductivity arises from 
strong correlation effects since the many-body energy lev- 
els favor electrons residing on the same molecule as op- 
posed to different molecules, resulting in a negative pair- 
binding energy. However, extensive numerical worki^ has 




FIG. 1: (Color online.) Molecular structures of C20 isomers: 
(a) cage, (b) bowl, and (c) ring. 



shown that the pair-binding energy in Cgo materials is 
likely positive for U/t < 4.5 with larger values of U/t 
inaccessible due to an increasingly severe sign-problem. 
Comparatively, C20 has a larger curvature and correla- 
tion effects measured in terms of U / 1 are therefore likely 
much more important since the curvature will decrease 
t and thereby increase U ft. We have previously stud- 
ied electron correlation effects in the C20 fullerene cageii 
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starting from a one-band Hubbard model denned in Eq. 
((T|). For the whole range of < U/t < oo, we find that 
the pair-binding energy is always positive ruling out the 
possibility of superconductivity induced solely by an elec- 
tronic mechanism. With increasing U/t there is, however, 
a very interesting metal- insulator transition predicted 11 
to occur around U/t ~ 4.1 in molecular solids formed 
with the C20 fullerene. Experiments have shown evi- 
dence for such solid phases of We expect that most 
of the strong correlations effects are fully captured by a 
one-band Hubbard model with a single uniform hopping 
integral. In order to develop predictive models of C20 
materials and to determine whether such materials will 
display metallic or insulating behavior it is then crucial 
to estimate how large U/t is for the C20 molecule, which 
is the main purpose of this paper. 

Mean field density functional theory (DFT) calcula- 
tions have yielded conflicting predictions of the relative 
stability of C20 isomers. This suggests that electron cor- 
relations, which are only approximately treated in the 
DFT calculations, could be very important in the pre- 
diction of the electronic structure of C20 isomers even 
though the geometry of the different isomers is fully cap- 
tured in the DFT calculations. We are, therefore, in- 
spired to study the isomers with an effective Hamiltonian 
approach, where a Huckcl Hamiltonian is complemented 
with an on-site electron correlation term, i.e., the one- 
band Hubbard model defined on a single C20 molecule 
as 

H = - >^ Ujic'^Cja + h.c.) + U / ]n^nii + £ p ^n;, 

(ij)tj i i 

(1) 

where c ia (ci a ) is an electron creation (annihilation) op- 
erator with spin a on site i, indices i,j run over 20 sites 
of the isomers, tij is the hopping integral between near- 
est neighbor (NN) carbon atoms i and j, U is the on-site 
electron correlation, ni = + is the number of elec- 
trons on site i, and e p is the on-site electronic energy due 
to the core ion and electrons in the carbon atom. Ac- 
cording to Ref. [H, we set e p = —8.97 eV. This on-site 
electronic energy term is a constant as long as the to- 
tal electron number in the molecules is fixed, and is not 
important in calculating, e.g., the pair- binding energy jii 
but needs to be taken into account when calculating the 
affinity energy of the molecules, defined as 

AE = E {20) - E (21), (2) 

where E(N) is the internal energy of the molecule filled 
with TV electrons from 2p atomic orbitals. 

The idea of fitting a tight-binding Huckel Hamiltonian 
to DFT energy levels for a fullerene molecule is not new 
and has been employed in, e.g., Ref. Il3lll4 However, 
inclusion of an on-site Coulomb interaction in the tight- 
binding model for the fullerene molecule has not been 
very well studied due to the difficulty of performing reli- 
able calculations in the presence of such a term. In par- 
ticular, the question of what value the on-site interaction 



U should take has not been answered for C20 isomers. 
In light of the metal-insulator transition predicted^ to 
occur around U/t = 4.1 a correct determination of U/t 
is crucial for modeling C20 based materials. Addressing 
this question is our main goal here. The paper is orga- 
nized as follows. 

First, since we want the tight-binding t term (Huckel 
Hamiltonian) to reflect the geometric character of each 
of the isomers, which is contained in the DFT energy 
levels, we fit the Huckel energy levels to energy levels 
obtained from DFT. This allows us to determine the 
effective value of t for the three isomers. Here we as- 
sume uniform hopping integrals t inside the cage and 
bowl molecules for the Huckel Hamiltonian, although 
within the DFT approach for the bowl isomer the hop- 
ping integrals are not uniform* even in the absence of 
any Jahn- Teller distortion. We expect this simplifica- 
tion for the Huckel approach to be of only minor impor- 
tance for the bowl and the cage. However, for the ring 
isomer, DFT calculations show that the ring structure 
is dimerized^ with alternating long and short bonds be- 
tween NN carbon atoms, which leads to a filled highest 
occupied molecular orbital (HOMO), i.e. an insulating 
molecule, in contrast to the uniform bonding case, where 
the HOMO is not completely filled (the molecule is metal- 
lic). We, therefore, consider two different hopping inte- 
grals for the ring Huckel Hamiltonian. As mentioned, 
inhomogeneous hopping integrals are always considered 
in the DFT calculations. We then study the effect of 
electron correlation in the neutral isomers by calculating 
the single-particle excitation spectra with different corre- 
lation strengths (U/t = 2, 3, 4, 5) and show the difference 
between quantum Monte Carlo (QMC) and DFT spec- 
tra for the neutral isomers. To complete our effective 
Hamiltonian approach, we proceed to estimate the on- 
site Coulomb interaction strength U by calculating the 
electron affinity energies of the isomers and comparing 
them to the experimental values. This allows us to de- 
termine U/t for each isomer and we can then study the 
effect of one-electron doping on the single-particle exci- 
tation spectra with QMC simulations. We compare the 
resulting QMC spectra to results obtained from DFT. 
These spectra should be directly comparable to possible 
future experimental PES and inverse PES spectra. 



II. ESTIMATION OF THE HUCKEL HOPPING 
INTEGRAL, t 

Different molecular geometries of C20 isomers de- 
termine different NN hopping integrals t in the Huckel 
description of molecules. In this section, we estimate 
approximate values of i's for the cage, bowl and ring 
structures of C20 molecules, by comparing tight-binding 
Huckel energy diagrams (U — 0) with the energy lev- 
els from DFT calculations. In DFT calculations for the 
cage and bowl, 2s, 2p x , 2p yi and 2p z atomic orbitals are 
considered, and the calculations are performed with the 
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TABLE I: Comparison of energy levels around the Fermi en- 
ergy from DFT calculation and Huckel Hamiltonian on the 
C20 isomers. The energy from DFT is in units of eV. D is 
degeneracy degree of the corresponding energy level, and N e 
is the number of electrons occupying these energy levels. For 
each isomer, the energy levels are listed in descending order. 

widely-used ABINIT codeJ^ Since there are 4 orbitals 
per carbon atom in the DFT calculation, the resulting 
energy levels are a mixture of a and ir bonds. The a 
bonds have either very low or very high energies. En- 
ergy levels around the Fermi energy mainly consist of ir 
bonds, which we use to construct the Huckel Hamilto- 
nian. The comparison of energy levels between DFT and 
tight-binding Huckel Hamiltonian should, therefore, be 
made around Fermi energy. Table |T] shows such a com- 
parison. To fit the Huckel hopping parameters t, we set 
equal the smallest energy gap of the tight- binding Huckel 
Hamiltonian with an equivalent gap in the DFT results. 
The degeneracy of the levels are here important and since 
the variation in the hopping integrals in the DFT ap- 
proach in some cases will split levels the equivalent gap in 
the DFT approach is not necessarily the smallest gap. In 
studies of strong correlation effects it is highly desirable 
to use the simplest possible model. Since for the fullerene 
cage the DFT and the tight-binding model with uniform 
hopping at U — both predict a half filled HOMO level 
and metallic behavior, we have not found it necessary to 
include variations in the hopping integral. The same is 
true for the bowl where a filled HOMO level is found. 
However, as previously mentioned, for the ring the two 
approaches disagree and we are therefore forced to in- 
clude a variation in the hopping integrals for the ring as 
we describe in detail below. 

We begin by discussing the fullerene cage. Here we find 
a half-filled 4-fold degenerate HOMO level and hence a 
metal in both the DFT (-0.98560 eV) and tight-binding 
(Oicago) results. In the tight-binding approach the first 
excited level is 4-fold degenerate at 2i cage . Assuming 
that this 4-fold level is split into a high lying 3-fold and a 
lower lying non-degenerate level due to slight variations 
in the effective hopping integrals in the DFT approach, 
the equivalent gap in the DFT approach should be from 
the 4-fold level at —0.98560 eV to some average of the 



3-fold and non-degenerate excited levels. In the extremal 
case we neglect the non-degenerate level and we then see 
that this gap should be close to the (1.72766 + 0.98560) 
cV= 2.71326 eV from which we infer that t cagc ~ 1.36 eV. 
Here we have taken the gap in the DFT calculations to go 
from the 4-fold degenerate level to the 3-fold degenerate 
level. Given the splitting of the excited levels in the DFT 
approach it seems plausible that this is a maximal value 
for the gap. However, in the absence of analytical results 
for the splitting of the levels in the DFT approach we note 
that this is not an exact bound. The estimate i cagc ~ 1.36 
eV is then likely an upper bound on the hopping integral. 
An extreme lower bound on t cage can be obtained if we 
use the smallest gap of (0.37579 + 0.98560) eV= 1.36139 
eV in the DFT approach, yielding ~ 0.68 eV. Given 
the large splitting of the non-degenerate level from the 
3-fold level this lower bound seems very unlikely to be 
attained and a more realistic value for i cagc is likely close 
to the upper bound of 1.36 eV which we mainly focus on 
in the following. 

We next turn to the bowl where one in both the DFT 
and tight binding approach (U — 0) finds a filled HOMO 
level and hence an insulator. For the bowl we com- 
pare the tight-binding gap of (0.47725 + 0.73764)t b owi = 
1.21489ibowl horn a 2- fold degenerate level to another 
2-fold degenerate level with the (-2.15325 + 5.02272) 
eV=2. 86947 eV gap in the DFT results yielding t howl ~ 
2.36 eV. As for the fullerene cage, we have here assumed 
that both 2-fold degenerate levels are slightly split in 
the DFT approach and we have used the largest rea- 
sonable value for the equivalent gap in the DFT ap- 
proach. Following the discussion of the hopping inte- 
gral for the fullerene cage we again expect the estimate 
ibowi — 2.36 eV to be an upper bound. In this case, 
a reasonable lower bound on the hopping integral can 
be obtained by taking the smallest gap in the DFT ap- 
proach of (-2.50728 + 5.00612) eV=2.49884 eV, resulting 
in tJJo™! — 2.06 eV, a relatively modest variation from our 
previous maximal estimate. 

Finally we turn to a discussion of the ring molecule. 
As mentioned, in this case a uniform tight binding model 
would predict a half-filled HOMO level where as the DFT 
approach shows a filled HOMO level. We are therefore 
forced to include a staggering in the hopping integrals 
in the tight-binding approach. For the ring DFT calcula- 
tion, we first generate a pseudo potential^ for the carbon 
atom, with 4 electrons (ls 2 2s 2 ) in the core state and 2 
electrons (2p 2 ) in the valence state. The resulting pseudo 
potential are then fed to the ABINIT^ to carry out DFT 
calculations on the 2 valence orbitals. On the tight- 
binding calculation side, we use 2 hopping integrals ti 
and t s to represent the hopping integrals of the alternat- 
ing long (bi = 2.609 Bohr)£ and short (b s = 2.260 Bohr)£ 
bonds, respectively. Let the average of the two bonds be 
b = (bi + b a )/2 = 2.435 Bohr. We then parametrize^ 
the two hopping integrals as t x /t r i ng = 1 — (b x — b)/b, 
where x — 1,8, and i r i ng is the average hopping inte- 
gral. With this parameterization we find ti = 0.928i r i ng 



and t s = 1.072f r j ng . From Table Q] we see that both 
DFT and tight-binding calculations now give an insu- 
lating molecule. If we again compare the tight-binding 
gap of (0.14400 + 0.14400)i ring = 0.288t ring from a non- 
degenerate level to another non-degenerate level with the 
gap between two non-degenerate levels in the DFT results 
of (-7.54251 + 8.45383) eV=0.91132 eV, resulting in an 
average t r i ng ~ 3.16 eV. 

We see that, compared to the other isomers, the hop- 
ping integral is significantly smaller in the fullerene cage 
isomer due to the large curvature of the molecule that 
reduces the NN overlap of the 2p orbitals. Consequently, 
i cag e for the C20 fullerene is also significantly smaller than 
what was found for the much bigger and less curved Ceo 
where one observes tc 6(l — 2.50 eV according to our cal- 
culation with ABINIT as well as Ref. [U or 2.72 eV ac- 
cording to Satpathy's early calculation^ It is also note- 
worthy that our estimate is likely an upper bound on 
^cagc- We also note that t, 2£ C age, which reflects the 
fact that i r i n g is the hopping integral of four 2p orbitals 
as apposed to £ C agc being the hopping integral of two 2p 1T 
orbitals. 



III. DOS FROM DFT CALCULATION 

Before we study the effect of electronic correlation 
with QMC, we calculate the density of states (DOS) for 
the neutral C20 isomers within DFT, shown in Fig. [2] 
We see that both the bowl and ring are insulators, with 
energy gaps = 2.4 eV and A°^t = 1.4 eV for the 

neutral molecule. The fullerene cage isomer is, however, 
metallic. Electron correlations are typically underesti- 
mated in DFT calculations and it therefore seems rea- 
sonable to assume that the inclusion of on-site electronic 
correlations would enlarge the gap, eventually turning 
the cage C20 isomer into an insulator. Such an effect was 
shown to occur in Ref. |Tl|. 



IV. EFFECT OF ON-SITE CORRELATION 

As mentioned, DFT calculations typically underesti- 
mate electron correlations. However, with the overlap 
integral, t, determined we can attempt to more closely 
describe the physics in the vicinity of the Fermi energy by 
using an effective one- band Hubbard model, Eq. (fTJ), to 
account for the electron correlations. To perform such a 
study of on-site electronic correlation on the C20 isomers 
we use the standard QMC algorithmj 17 ' 18 i 19 In regimes 
where the sign problem in this approach renders results 
unobtainable we have supplemented these results by ex- 
act diagonalization (ED) results. This allows us to deter- 
mine the DOS as a functions of U/t. The DOS is calcu- 
lated for each of the neutral isomers for U/t = 2, 3, 4, 5. 
Our results are shown in Fig. (3) As expected, the energy 
gap increases with increasing U/t values for the bowl and 
the ring. The dependence on U/t is clearly non-trivial. 
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0) (eV) 

FIG. 2: (Color online.) DOS from DFT for neutral cage, 
bowl and ring C20 isomers. Fermi energy is located at u> — 0. 
Shaded areas are occupied by electrons. 



For the fullerene cage the initially metallic molecule be- 
comes insulating with increasing (U/t) cage . For molec- 
ular solids formed out of this isomer a metal-insulator 
transition is therefore expected around (U/t) cage = 4.1»ii 
In the next section we estimate the on-site Coulomb in- 
teraction U using the affinity energy of the isomers. How- 
ever, assuming that the on-site electronic energy scale is 
U ^ 10 eV for all the isomers ; 20 ! 21 ' 22 ? 23 we can do a rough 
estimate of U/t and hence estimate the gap for the neu- 
tral molecule from the above QMC results for the DOS. 
Using this value for U we find that U/t ~ 4.2,3.2 for 
the bowl and ring isomers respectively using t r i ng ~ 3.16 
eV and the upper bound ibowi — 2.36 eV. For these 
values of U jt and from Fig. [31 we estimate energy 
gaps for the neutral bowl and ring isomers to be about 
Abowi = 2.0t b owi = 4.72 eV and A ring = 0.9i ring = 2.73 
eV, respectively. For the cage isomer (U/t) cage should 
be around 7 using the upper bound t cagc = 1.36 eV. Un- 
fortunately, the sign problem prevents us from calculat- 
ing the DOS and hence the gap by QMC simulations 
for (U/t) cagf} > 3 and we have to resort to the much 
more time consuming ED approach. Previously, using 
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FIG. 3: (Color online.) Evolution of the DOS with differ- 
ent U/t values for neutral C20 isomers: cage, bowl and ring. 
U/t — 5 for cage isomer is from exact diagonalization.— Oth- 
ers are from QMC simulations. Fermi energies are located at 
u = 0. 



ED techniques, the DOS and gap at (U/t) c 



5 has 



been determined— yielding a gap of 1.89 eV. (See Fig.[3j) 
We have repeated this calculation at (U/t) cagc = 10 
finding a gap of 7.67 eV. Linearly interpolating between 
these values we estimate the gap for the fullerene cage at 
(U/t) 



cage 



7 to be Acagc — 3.1£ ca g C 



4.2 eV. 



V. ESTIMATE OF THE ON-SITE 
INTERACTION, U/t 

Electron correlations are important for an accurate 
calculation of the affinity energy, as shown by early QMC 
simulations. See, e.g., Ref. (2J. Therefore, starting 
from the affinity energy, one can use QMC simulations 
to determine electron correlations, which in the case of 
Hubbard model are represented by the on-site interac- 
tion U. Since the electron affinity energies, AE, were 
measured experimentally,— we can use these energies for 
a more refined estimate of the value of the on-site ef- 




FIG. 4: (Color online.) Variation of affinity energy with U/t 
for a C20 cage, bowl and ring, from QMC and ED calculations. 
The experimental AE values are shown by the horizontal lines 
intercepting the AE ~ U/t curves. The corresponding U/t 
values are shown by the vertical lines intercepting the curves. 
Results are shown using the upper bound for t for the cage 
and bowl. 



fective Coulomb interaction U for each of the isomers. 
Ground state energies of the isomers are calculated with 
projection QMC technique^ in the Hilbert space of fixed 
particle number N and magnetization and converted to 
electron volts using the previously determined estimates 
for t. The affinity energy AE = E(20) - E(21) is then 
subtracted and shown in Fig. 2] as a function of U /t. 
Results are shown using the upper bound for t for the 
cage and bowl. The sign problem prevents us from sim- 
ulating larger U/t values with QMC. Instead we use ED 
to obtain a series of AE for U/t = 2, 5, 6, 7, 8, 10 for the 
cage isomer. Note that, for the fullerene cage the elec- 
tron affinity changes slope around the expected critical 
value of (U/t) cagc — 4.1 «ii We see that all 3 curves are 
approximately linear around the experimental affinity en- 
ergy values^ (AE cage = 2.25 eV, AE how \ — 2.17 eV, and 
2.44 eV). Thus, a linear interpolation gives 



AE Y [ n g 

t^cagc — 7.1t ca g C , t/bowl — 4.30£bowl 



and U^'mg — 3.27£rin 



If we for the fullerene cage perform a similar analysis us- 
ing i™ e = 0.68 eV we find instead an even larger value 
for (£//t) ca ge — 12.2. Likewise, we find for the bowl using 
= 2.06 eV a value of (U/t) howl = 4.9. 
The above results (Fig. [4j) show that for cage, bowl and 
ring isomers the on-site Coulomb interaction energies are 
U cage = 7.1i cagc = 9.67 eV, U how i = 4.30t bo wi = 10.15 
eV, and U T i ag = 3.27i r i ng = 10.33 eV, all of which are 
reasonably close to 10 eV, i.e., the value that was used 
in the previous section. The variation of U ft between 
the isomers is minor and the consistency between these 
results and those of the previous section supports our 
approach of determining t from DFT calculations. Using 
the upper bounds on U/t for the cage and the bowl with 
the associated minimal values for the hopping integrals 
results in U = 8.3 eV and U — 10.09 eV, respectively. It 
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is noteworthy that the value of (U/i) ca , ge is clearly beyond 
the metal-insulator transition point U/t = 4.1 predicted 
in Ref. [ll|. This suggests that undoped molecular solids 
formed of dodecahedral C20 are insulators, and the iso- 
lated molecule is likely not Jahn- Teller active^ 

In conclusion, the estimation of U/t using the molecu- 
lar affinity energy in combination with a determination of 
the hopping integral t from DFT is a reasonable approach 
for the cage, bowl and ring isomers. We summarize all 
the isomer parameters in Tab. [TTJ 



VI. ENERGY GAPS AND EFFECT OF 
ELECTRON DOPING 

We proceed to study the effect of one electron dop- 
ing on the single particle excitation spectra of C20 bowl 
and ring isomers with QMC simulations, using the hop- 
ping integrals tbowi — 2.36 eV and i ring = 3.16 eV 
and on-site Coulomb interaction strength U estimated 
above. A similar study for the fullerene cage would re- 
quire time consuming ED calculations which we have not 
performed. The DOS are shown in Fig. [5] We see 
that with one electron doping, spectral peaks for both 
bowl and ring move toward Fermi levels. Clearly, the 
one-electron-doped bowl becomes metallic. For the one- 
electron-doped ring, the spectral weight is very close to 
the Fermi energy. However, within the precision of the 
QMC results, we expect the doped ring to remain an in- 
sulator with a gap A^ MC = 3.2 eV. This is in contrast 

ring 

to DFT energy levels, Table U which show that with one 
electron doping both the bowl and the ring isomers be- 
come metallic. 

From these figures we can also estimate the energy gaps 
for the neutral bowl and ring to be around A? M 9 = 

t -* DOW 1 

4.0 eV and A^g = 3.6 eV in reasonable agreement 
with the rough estimate given in section HVl We note 
that these values are much larger than the corresponding 
DFT values A£ F J = 2.4 eV and A^ g T = 1.4 eV. From 
the above discussion we expect our previous estimate of 
the gap for the fullerene cage of A cagc = 3.1£ cagc = 4.2 
eV at (U/t) ca%c ~ 7 to be relatively precise. Hopefully, 
these estimates of the gaps can be compared to future 
experimental PES and inverse PES of the neutral C20 
isomers. 





t (cV) 


U (eV) 


U/t 


cage 


0.68-1.36 


8.3-9.67 


7.1-12.2 


bowl 


2.06-2.36 


10.09-10.15 


4.30-4.9 


ring 


3.16 


10.33 


3.27 



TABLE II: Hopping t and on-site Coulomb interaction U for 
C20 isomers: cage, bowl and ring. 




to(eV) 

FIG. 5: (Color online.) QMC results for the DOS of one- 
electron doped C20 (solid lines) compared with the neutral 
molecules (dashed lines). For the bowl, (U/t)bowl = 4.30; 
for the ring, (J7/t) r ing = 3.27. Shaded areas are occupied 
by electrons. Fermi energies are located at lo = 0. Energy 
units have been converted to eV using ibowi = 2.36 eV and 
iring = 3.16 eV, respectively. Energy gaps for the neutral 
molecules and doped ring are shown in the figure. 



VII. CONCLUSIONS 

We proposed an effective Hamiltonian approach to 
study the electronic correlations in C20 isomers with 
QMC simulations and exact diagonalization. The hop- 
ping integral, t, in an effective one band Hubbard model 
Hamiltonian are determined by comparing DFT energy 
levels with the tight-binding Huckel energy levels. On- 
site Coulomb interactions, U, are then determined by 
comparison to the experimental affinity energies of the 
isomers. With these estimated parameters, QMC calcu- 
lations of the resulting effective Hubbard model then pre- 
dicts insulating behavior of the neutral (cage, bowl and 
ring) and one-electron-doped (ring) isomers and metal- 
lic behavior of the one-electron-doped bowl isomer. We 
find qualitative agreement between QMC and DFT cal- 
culations for neutral (bowl and ring) isomers and one- 
electron-doped bowl isomer, although QMC gives much 
larger energy gaps for the neutral isomers. For the neu- 
tral cage isomer, the QMC prediction (insulating) is qual- 
itatively different from DFT calculations (metallic), since 
the cage C20 is the most strongly-correlated molecule of 
the three isomers with a ratio of on-site Coulomb inter- 
action and hopping integral of (U/t) cagc ~ 7.1 — 12.2 ex- 
ceeding the value of 4.1 for the predicted metal-insulator 
transition. Results presented in the paper await a com- 
parison with possible future PES and inverse PES exper- 
iments on the gas phase of C20 isomers. 
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